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Abstract 

Dispersive corrections to the total cross section for high-energy scattering from a heavy nucleus 
are calculated using a matrix model, based on the triple-Pomeron behavior of diffractive scattering 
■ from a single nucleon, for the cross section operator connecting different states of the projectile 

m 

nucleon . Energy-dependent effects due to the decrease in longitudinal momentum transfers and 
the opening of more channels with increasing energy are included. The three leading terms in an 
expansion in the number of inelastic transitions are evaluated and compared to exact results for 
the model in the uniform nuclear density approximation for the the scattering of nucleons from 
O , Pb 20S for laboratory momenta ranging from 50 to 200 GeV/c . 
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I. INTRODUCTION 



One consequence of the composite nature of nucleons is a decrease in nucleon-nucleus 
total cross sections due to transitions between different internal states of the projectile 
nucleon. This decrease can be calculated 0] using an operator to represent a generalized 
nucleon-nucleon total cross section, with the matrix elements representing the probability 
amplitudes for forward scattering transitions between different states of the nucleon. 

The general problem of the propagation of composite particles through a nucleus has a 
long history || f|, [^, || [7|, R] . If the longitudinal momentum transfer due to the different 
masses of the different states of the composite system cannot be ignored the calculation of 
the dispersive corrections becomes quite complicated. Because of the difficulty of this calcu- 
lation, and uncertainties in the nature of the cross section operator, only the leading term in 
an expansion in the number of inelastic transitions has been evaluated §J . (A closely related 
effect in quasi-elastic electron scattering can complicate the analysis of color transparency 
effects in these reactions ||10|| .) In this paper we test the accuracy of this approximation, and 
the convergence of the expansion, using a simple finite matrix model for the cross section 
operator and taking the uniform nuclear density limit for which the exact result, including 
all orders in transitions, can also be calculated. This matrix model is consistent with the 
known triple- Pomeron behavior of high-energy, high-mass single diffraction dissociation from 
a single nucleon, but otherwise has a high degree of arbitrariness: it is used only for lack of 
a reliable microscopic model for the internal degrees of freedom of highly excited nucleons. 
It has some features in common with a "simplified example" used by Van Hove |7|] for the 
limit of zero longitudnal momentum transfer. 

Section II reviews the transition expansion for the dispersive corrections to the nucleon- 
nucleus total cross section, including the effects of the longitudinal momentum transfers due 
to the different masses of the the different excited states of the nucleon. In Sect. Ill it is 
shown that the terms in this expansion simplify considerably in the uniform nuclear density 
limit, with each term in the expansion represented as a sum over products of transition 
amplitudes weighted by a function of the differences among cross sections and differences in 
the masses of the different nucleonic states. In the uniform density limit it is also possible 
to write the exact result, including all orders in the number of transitions, in terms of the 
exponential of a single position-independent operator, as shown in Sect. IV. The represen- 
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tation of the transition operator by a finite matrix, with a dimension which increases with 
energy, is discussed in Sect. V. This matrix is chosen to have a form consistent with the 
triple-Pomeron behavior of the nucleon-nucleon single diffraction dissociation, but is other- 
wise highly arbitrary. Using this matrix the formulas are evaluated and the results presented 
in Sect. VI. The results are summarized and discussed in the concluding Sect. VII. 



II. EXPANSION IN INELASTIC TRANSITIONS 



It has long been known that a version of the eikonal approximation holds for an infinite 
but restricted class of Feynman diagrams which includes inelastic transitions among different 
states of the projectile |TT| . The result is equivalent to the eikonal approximation in coupled- 



channel potential theory [p^ , 13], |T4|j , and leads to an expression for projectile-nucleus total 
cross sections which can be written as 

a{A) Total = 2Re J d 2 b(A\ (l\f (b, {r a })\l)\A), (1) 

where \A) is the ground state of the nucleus and |1) is the lowest mass eigenstate of the 
projectile system. The profile function f is an operator in the internal space of the projectile: 

f(b, {r Q }) = 1 - Zexpi-i £ f dzu a (v)], (2) 

where Z indicates a z-ordered product and 

u a {v) = {m/ Pl )e-^v{v-v a )e^\ (3) 

with v(r — r a ) the effective potential operator produced by a static target nucleon at r a . 
Here p is the longitudinal momentum operator, diagonal in the mass eigenstates of the 
projectile, with 

p j ^p 1 -(m^-m 2 1 )/{2p 1 ), (4) 

where pi is the initial momentum of the projectile in the laboratory system where the nucleus 
is at rest and rrij is the mass of the jth excited state of the projectile. 

Assuming the different v a do not overlap, ignoring nuclear correlations, and assuming a 
large nucleon number A, the expression for T simplifies to 

f(6) = (A\T(b); {r a }\A} « 1 - Zexp[-{A/2) J dz lP (b, z 1 )a{z 1 )\, (5) 



where p is the nuclear density, normalized to one, and 

<j{z\) = exp{— ipzi)aexp(ipzi), (6) 

with a the cross section operator M for scattering of the projectile from a single target 
nucleon. Assuming that the corresponding elastic and diffractive scattering amplitudes are 
purely imaginary, the total projectile-nucleon cross section is 

a T = (l\a\l), (7) 

while the cross section for single diffraction dissociation of the projectile interacting with a 
single nucleon, summed over all diffractively excited states of the projectile, at momentum 
transfer squared t — 0, is 

da diff /dt = [(l|a 2 |l) - (l|a|l) 2 ]/(i67r). (8) 

If 

r(6) = (l|f(6)|l), (9) 

then 

a{A) Total = 2Re J d 2 bT{b). (10) 

Since the dispersive corrections are small compared to the total cross section it is useful 
to separate Y{b) into two parts: 

r(6) = r G (6)-r D (6). (11) 

Defining the dimensionless absorption parameter 

t{b) = {A/2)a T J dzp{b,z), (12) 

the main Glauber contribution, which does not include dispersive corrections, is simply 

T G (b) = l-exp[-t(b)], (13) 

while the diffractive correction is given by 

T D (b) = (l\Zexp[-(A/2) J dzp(b,z)a(z)]\l) - exp[-t(b)]. (14) 

Previous calculations taking into account longitudnal momentum transfers of 
have included only the leading second order term in an expansion in the number of inelastic 



transition. While this is almost certainly accurate for light nuclei, it is not clear whether or 
not it is adequate for heavy nuclei. Below we develop expressions for the general terms in 
the expansion and evaluate them in a simple but possibly not completely unrealistic model. 

We begin by separating the cross section operator a into its diagonal and off-diagonal 
(transition) parts: 

a = cr d + a t , (15) 

where 

(i\*d\j) = (16) 

so that a t has only off-diagonal matrix elements. We now expand in powers of the 
transition operator: 

oo 

M6) = Er& } (6)- (17) 

n=2 

Here the leading n = 2 term contains contributions from processes in which the nucleon 
makes two transitions: one from the ground state to a higher-mass state, then another back 
to the ground state. The next (n = 3 ) term contains the contributions from processes with 
three transitions, with the two intermediate projectile states being neither the ground state 
nor equal to each other. 

Using a derivation analogous to that for time-dependent perturbation theory |TJ|] one can 
show that the F£ \b) are given by the z-ordered integrals 

, s /-OO fZ2 



r^(6) = (-A/2) n e / dz n p(b, z n )... dz lP (b } z x ) 

J —oo J —oo 

(l\a td (z n ) . . . o- M (2;2)^ d (zi)|l), (18) 

a td (z) =U(z)- 1 a t U(z), (19) 

U(z) = exp[-(A/2) [ Z dz xP {b,z x )o d + iyz\ (20) 
J —00 

is diagonal in mass eigenstates. These expressions can be simplified if we replace z by the 
dimensionless variable 

u(b,z) = (A/2)a T [ Z dz lP (b,z 1 )/t(b), (21) 

J —00 



with 



where 
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so that u(b, — oo) = and u(b, oo) = 1. If we also define the dimensionless cross section 
operators 

x = a/a T , (22) 
x d = a d /a T , (23) 
x t = a t /a T , (24) 

and, for future use, 

x = x-l. (25) 

We can then write 

T%\b) = {-t{b)) n exp[-t{b)\ f 1 du n ... r du 2 r d Ul 

Jo Jo Jo 

(l\x t d(u n )...x t d(ui)\l), (26) 

where 

x td (u) = U(z(b, u))- 1 x t U(z(b, «)), (27) 

with z(b, u) the inverse of u(b, z) for fixed b. 

Inserting complete sets of mass eigenstates between the x td (u) operators in Eqn. (20) 
gives 

r& } (&) = [{-t{b)) n /n>]exp[-t{b)\ £ (i|^bn-i)...0 2 |^b'i)0'i|^|i) 

31,— Jn-l 

/ {n) (^Jl,J2,...,Jn-l), (28) 



where 



r n) (b;j 1 ,j 2 ,---,jn-i) = n\ du n ... / du 2 / dui 

JO JO JO 

exp[t(&)u n (xi - x jn _J - iznfa - p jn _J] 
. . . exp[f(&)it 2 (:Cj a - XjJ - iz 2 (p j2 - p h )\ 
exp[t(b)u 1 (x h - xi) - iz l (p h - p^], (29) 

with Zj = z(b,Uj) and Xj =< j\x\j >. Since x t has only off-diagonal matrix elements, terms 
with equal successive ji's do not contribute to the sum in Eqn 22. 

Clearly the functions defined above depend only on the differences between successive 
x'jS and p'jS, and thus are unchanged if these variables are replaced by 

Xj =Xj-\ (30) 



and 

Pj=Pj-Pi- (31) 

Then can be written 

f (n) (b;j 1 ,j 2 ,---,jn~i)=n\ du n ... du 2 du x 

Jo Jo Jo 

exp\-tx jn _ x {u n - Un-i) + ivj n _ x {z n - z n _i)] 

. . . exp[-tx jl (u 2 - ui) + ip j2 (z 2 - z-l)], (32) 

where, because of the u (and z) ordering, the differences between u's and z's in the paren- 
thesis are never negative, and the real parts of the 's, needed to calculate the total cross 
sections, are always less than their limits as the pj 's approach zero. Each exponential in 
this expression acts as a propagator for the projectile from the location of one transition 
to that of the next and includes absorptive and phase-changing parts depending upon the 
state of the projectile at this stage of its journey through the nucleus. 

III. UNIFORM DENSITY LIMIT 

The evaluation of the expressions for the 's is complicated by the fact that z and u 
are in general not simple functions of one another: the relation between them is determined 
by the shape of the nuclear density function p and depends upon the impact parameter. 
For heavy nuclei p is well-approximated by the simple Woods-Saxon form, but the relation 
between z and u is still not simple. For the heaviest nuclei the surface thickness is much 
less than the nuclear radius so that it may not introduce excessive errors to replace p by 
its uniform density limit, especially at small momentum transfers and , in particular, for 
evaluating the dispersive contribution to the total cross section. In this limit 



giving 



where 



p(r)=p Q(R-r), 

t(b) = Aa T p V R 2 — b 2 Q(R - b), 
Po = l/(4rf/3) 



(33) 

(34) 
(35) 
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and R ~ r^A 1 ^ , with r ~ 1.14/m . (With this expression for t{b) there is an analytic 
expression [|j] for the main Glauber contribution to o(A)Totai-) In this limit z and u are 
linearly related: 

z = VR 2 -b 2 {2u- 1) (36) 
for < u < 1. The arguments of the exponentials in Eqn. 26 then simplify considerably: 

- ZjJ(b)(u n +i ~ u n ) + ip jn (z n+1 - z n ) = -y jn VR 2 - b 2 (u n+1 - u n ), (37) 
where the complex number 

yj = XjAa T p - 2ip j (38) 

is independent of b and the u's . In the uniform density limit, then, the effects of longitudinal 
momentum transfers are taken into account simply by adding an imaginary part to each 
diagonal matrix element Xj, and modifying the calculations of the functions accordingly. 

IV. EXACT RESULT 

In the uniform density limit an exact expression for the dispersive correction correction to 
the total cross section can be found in terms of the exponential of a z-independent operator. 
The simplest derivation of this result starts with Eqn. 5 and removes the z-dependence of 
a(z) by adding a term proportional to p to the operator in the exponent. For a given impact 
parameter a single z-independent matrix is involved so that the z-ordering in Eqn. (2) can 
be ignored, giving 

r D (6) = (l\exp(-M(b)\l) - exp(-t(b)), (39) 
with the z-independent operator 

M{b) = [Ap a - 2t(p - Pl )]VR 2 -b 2 . (40) 

This can also be written as 

F D (b) = exp(-t(b))[(l\exp(-M(b))\l) - 1], (41) 

where 

M(b) = M{b) - t(b)l = [A Po a T x - 2tp]VR 2 -b 2 . (42) 
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Either of the expressions (33) or (35) can be evaluated by expanding the exponential of the z- 
independent operator in a power series , with the second converging somewhat more rapidly. 
(Using the matrix model below one has to include of the order of 50 terms in the expansion, 
and there is considerable cancellation, so the individual terms must be calculated to high 
accuracy.) 

Another approach for evaluating To(b) depends upon the fact that, since the nucleon 
always enters the reaction in its ground state, the full operator exp(-M) is not needed: it 
is sufficient to work with the reaction-modified state 

\V) = exp{-M)\l), (43) 

with 

T D (b)=exp(-t(b))[(l\V)-l]. (44) 

Expanding the exponential 

oo 

\V) = Y,\Vn), (45) 

n=0 

where 

\Vn) = (l/d)(-M) n |l) = -(l/n)M\V(n - 1)) (46) 
can be calculated recursively starting with |V0) = |1). 



V. MATRIX MODEL 

The expressions for above involve the matrix elements of the dimensionless cross sec- 
tion operator x and the longitudinal momentum operator p . In this section a model for these 
operators is developed which, although highly arbitrary, is consistent with the experimental 
high-energy behavior of diffraction dissociation, which is in turn approximately consistent 
with that expected from the leading triple-Pomeron behavior. Ignoring contributions from 
secondary Regge poles, and taking the Pomeron intercept a(0) = 1, this leads to the simple 



behavior at momentum transfer squared t — |16| , |17 

da/(dtdM 2 ) = a T 3/ Vpp/(16vrM 2 ), (47) 



where <jt is the nucleon- nucleon total cross section and gppp ~ 0.364m6 1 / 2 @ is the 
triple-Pomeron vertex, while the mass-square of the diffractively excited nucleon runs from 
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Mmin ~ 1.5(GeV/c 2 ) 2 to M^ ax m m\ + 2pim 7r , the later condition following from the 
requirement that the longitudinal momentum transfer be less than m^, the inverse of the 
range of the strong force. (At very high energies the effective Pomeron intercept is greater 
than one and the above must be modified. Details can be found in JB5], 



The continuous range of m 2 between M^ in and M^ ax can be approximately replaced 
by a finite number of "bins" of width Am 2 centered at m 2 . The operator x can then be 
represented by a finite matrix with elements (j\x\i) constrained by 

daj/dt « 4(j|x|l) 2 /(167r) « [dcr / (dtdM 2 )] m 2 Am 2 , (48) 

or 

(j|x|l) 2 « g PPP Am 2 /(a)( 2 m 2 ) (49) 

To complete the model one must also have a prescription for the sizes of the mass bins. For 
simplicity equal spacing in m 2 is used below: 

m) = m\ + rr%{j-l), (50) 

with m a parameter determined by the spacing of the low energy diffractively produced 
resonances. Eqn. (35) then takes the simple form 

(j\x\l) 2 « (gppp/o-T /2 )/[{mi/m ) 2 +j - 1). (51) 

With the expression above for M^ ax , the dimension N of the matrix is given by 

N w 2p 1 m 7T /m 2 (52) 

which increases linearly with pi , the momentum of the incident proton in the rest frame of the 
target nucleus. This expression for N is clearly only a rough estimate, but fortunately excited 
states with j near iV do not contribute much to the dispersive correction compared to lower 
states: changing N slightly does not affect the results below appreciably. Unfortunately 
diffraction dissociation constrains only one row (and column, from the assumed symmetry) 
of the matrix (j\x\i) . For small m? and large m 2 one can argue that the triple- Pomeron 
behavior should still be valid and the Am 2 /m 2 dependence should still hold. For simplicity, 
here it is assumed that every off-diagonal element of the matrix is given by 

(Mi) = \I(9PPp/[4 /2 (a 2 + \j-i\)}, (53) 
10 



where a = mi/m . This expression is consistent with both experiment and the triple- 
Pomeron behavior for i — 1 and large j, but is only a guess elsewhere, especially when i and 
j are comparable in size. Furthermore, although the matrix is in general complex, it will 
below be assumed real. This is done mainly because the phases of the matrix elements are 
unknown (except for that of rr n ) and is consistent with the fact that the real part of the 
forward proton-proton scattering amplitude at high energies is known to be small. 
Finally, we need an expression for the diagonal elements Xj = (j\x\j) : 

Xj = l + d(j-l), (54) 

which allows the cross section for nucleons to scatter from excited nucleons to increase 
with the degree of excitation. This means that highly excited states will be absorbed more 
strongly than lower states while propagating between transitions. (If d — all diagonal 
elements of x are unity and and the expressions above simplify considerably.) It would 
probably be more reasonable for Xj to approach a limiting value as j increases, but this 
would introduce still more parameters and assumptions into the model. 

With these assumptions both Xj and are linear in j, and so are the complex numbers 
Uj defined in Eqn. 34: 

Vj = U - l )V2, ( 55 ) 

where 

y 2 = Aa T p d + iml/px. (56) 

For small j's the influence of the longitudinal momentum transfer, given by the imaginary 
part of z/2, decreases as 1/pi, but the dimension of the matrix increases as pi so that for the 
heaviest excited nucleon included we have 

y N = 2Aa T p dm 7V p 1 /ml + i2m n , (57) 

with an imaginary part which is independent of p\ and a real part which increases linearly 
with momentum. 

VI. RESULTS 

The formulas above have been evaluated for scattering of nucleons from 208 Pb for incident 
laboratory momenta from 50 to 200 GeV/c, and for n, the number of inelastic diffractive 
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TABLE I: Dispersive reductions in the total cross sections for scattering of high-energy nucleons 
from Pb 20S . The figures in parentheses are obtained if the longitudinal momentum transfers are 
set to zero. 



Lab 
Mom 
GeV 


Dispersive Reductions in <JTotal{A = 208) (mb) (ApL = 0) 


Order in Inelastic Transitions 


Sum 


Exact 


2 


3 


4 


50 


45.4 (57.9) 


-16.4 (-39.2) 


5.6 (31.4) 


34.6 (50.1) 


34.2 (37.1) 


100 


59.6 (66.9) 


-40.6 (-61.4) 


29.6 (65.7) 


48.6 (71.2) 


37.6 (38.5) 


150 


66.0 (70.6) 


-56.9 (-73.0) 


54.4 (88.2) 


63.5 (85.8) 


38.5 (38.9) 


200 


69.7 (72.9) 


-68.2 (-80.6) 


75.3 (104.8) 


76.8 (97.1) 


38.9 (39.2) 



transitions, ranging from 2 to 4. The corresponding reductions in the total cross sections 
can be written as 

aff = J db4nbReT ( v\b). (58) 

Using gppp = 0.363mb 1 ^ 2 and = 38.5m6 gives the coefficient in Eqn. 39 (gppp/cr^' 2 ) 1 ^ 2 = 
0.24. We also take tUq = 1.5GeV 2 , a 2 = 0.5 and d = 0.1, although any values of the same 
order of magnitude would be just as reasonable. 

The shapes of the integrands for some 50 and 200 GeV are shown in Figs. 1 and 2, 
respectively, while the values of the diffractive reductions in the total cross sections are given 
in Table I. The most surprising feature of these results is that, although the contributions 
are all small compared to the un-corrected Glauber cross section of about 2580m6 , at high 
energies the expansion in the number of inelastic transitions does not converge well at all. As 
noted above, as the energy increases more and heavier excited nucleons are included, and the 
correction at each order increases. Because of the poor convergence, however, the expansion 
in the number of inelastic transitions is not very useful, especially at higher energies: it gives 
only an order-of-magnitude for the exact result, with the individual terms oscillating in sign. 
The leading second order term in particular is always of the right order of magnitude but 
larger than the exact result, with the error increasing from 33% at 50 GeV/c to nearly 80% 
at 200 GeV/c. 

The influence of the longitudinal momentum transfer was studied by comparing the result 
calculated from the formulas above with those with the longitudinal momentum transfers 
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1 2 3 4 5 6 7 



b(fm) 



FIG. 1: Integrands for diffractive reductions in nucleon-P6 total cross sections as a function of 
impact parameter b at 50 GeV in 2nd through 4th order in the number of inelastic transitions. 
The dashed curve shows the sum of these corrections, while the heavy solid curve shows the exact 
complete correction. 

dropped (so that y 2 becomes a real number). The results in this limit are given as the 
numbers in parenthesis in Table 1. Although the longitudinal momentum transfers reduce 
significantly the magnitudes of the individual terms in the expansion, , they have a relatively 
small effect, decreasing with increasing energy, on the exact results. 
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FIG. 2: Integrands for diffractive reductions in nucleon-P6 total cross sections as a function of 
impact parameter b at 200 GeV in 2nd through 4th order in the number of inelastic transitions. 
The dashed curve shows the sum of these corrections, while the heavy solid curve shows the exact 
complete correction. 

VII. CONCLUSION 

Previous calculations of the dispersive corrections considered here have considered only 
a small number of channels, ignored the longitudinal momentum transfer, or included only 
the lowest order term in the the transition expansion. In the context of our model it has 
been shown that all of these can lead to large errors. It should be remembered, however, 
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that the matrix model used above has many arbitrary features, even though it is roughly 
consistent with what is known about about high-mass single diffraction dissociation from a 
single nucleon. In particular, the amplitudes for transitions from one highly excited nucleon 
state to another are essentially unknown, and the expressions used in the model are simply 
guesses based on the known behavior of the amplitudes for excitation from the nucleon itself. 
It would be useful to repeat the calculations with other assumptions for these amplitudes in 
order to get some idea of the dependence of the results on the assumptions. 

In the model used here the diffractive corrections to the total cross sections are all small 
compared to the total cross section itself, but their expansion in the number of inelastic 
transitions does not converge well at higher energies. In particular, the leading second order 
correction, which has been used to estimate the size of the diffractive correction, is too large 
by nearly 80% at a laboratory momentum of 200 Gev/c. 

It would be interesting to extend these calculations to single diffraction dissociation from 
nuclei, since for these processes there is no zeroth order term, corresponding to the large 
Glauber contribution to the total cross section , so that the corrections due to higher order 
terms might be relatively quite large. A preliminary investigation suggests that it should also 
be possible to do an exact calculation in this case in the uniform density limit. (A calculation 
of single diffraction dissociation from the deuteron would also be very interesting, and might 
put additional constraints on the assumptions that go into the matrix model.) "Coherent" 
diffraction dissociation, where the nucleus remains in its ground state, would be particularly 
simple to calculate, but probably experimentally challenging. One could also calculate a 
"nuclear inclusive" cross section in which all nuclear excited states are summed over. 
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